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1. Introduction 

Next yeai" CERN's Large Hadron Collider (LHC) will start delivering proton beams for physics 
collision. The LHCb experiment is designed to exploit the enormous LHC potential in the b-quark 
sector for measurements of the CKM parameters to such a high precision that possible contributions 
from TeV-scale New Physics to the mixing mechanism will become visible. To give an idea, the 
production of b hadrons at LHCb is expected with the annual yield of 10^^ b-b pairs [|]]. Possible 
future super-B factories would further extend the set of high precision b-physics measurements 

This programme can provide a stringent test of the Standard Model and potentially lead to 
the discovery of New-Physics only if at the same time a significant progress on the theory side is 
made. To get a flavor about the required precision it is useful to have a look at the experimental and 
theoretical situation for a few low-energy flavor-violating observables where non-Standard effects 
were expected to contribute. 

Let us start with the inclusive radiative B-meson decay. The world average performed by the 
Heavy Flavor Averaging Group [^] for Ey> 1.6 GeV yields the branching ratio 

3g{B^X,y) = (3.55 ±0.24^°:?^ ±0.03) x 10"^ , (1.1) 

to be compared with the Standard Model NNLO analysis of Ref. [Q], which for the same cut on the 
photon energy gives 

^(B^Xs7) = (3.15±0.23) X 10"^ . (L2) 

The values are consistent basically within one (combined) sigma (about 10%), which implies that 
the difference between the Standard Model (SM) and the experimental numbers can be of order 
20%. Notice that this estimate does not depend on theoretical inputs from the lattice. 

New Physics in principle can be found also in purely leptonic (or D^) decays, which can 
be enhanced by charged Higgs exchange contributions in any model with two Higgs doublets 
Again, the average of the experimental numbers for SS{B — > Tv) from Belle and Babar ||6[ ^ is in 
good agreement with the SM theoretical computation although the total error is above 30% in the 
first case and around 20% in the second. On the theory side this is mainly due to the uncertainty on 
Kib and on the decay constant F^, which will be discussed in more detail in the next section. The 
experimental error on the other hand is expected to decrease in the future. At a Super-B factory 
with hundred times the luminosity of Belle the branching ratio would be measured with a precision 
of 3%. 

The leptonic decays of the neutral Bj^) meson are very rare (the SM branching ratio is 0(10^^)) 
and they haven't been observed so far. The most recent experimental upper bound on 3S{B^ — > 
fj.^ fj.^) is 1 X 10^^ from CDF This decay is included among the LHCb physics goals with 
an SM expectation of 20 events per year There is quite some excitement around this channel 
as it can be significantly enhanced in various extensions of the Standard Model. For example 
in the Minimal Supersymmetric Standard Model with large tanjS (where tanjS is the ratio of the 
two neutral Higgs field vacuum expectation values) the enhancement can be up to three orders of 
magnitude compared to the SM. That is due to the appearance of flavor changing couplings of the 
neutral Higgs bosons generated by non-holomorphic terms after supersymmetry breaking [^. 

Finally New Physics might contribute to D — D mixing, which has been recently observed by 
Babar [|I0|]. It is hard to quantify the possible size of non-Standard effects here as the SM theoretical 



2 



SM parameters and heavy quarks on the lattice 



Michele Delia Morte 



predictions are very uncertain. It is also not clear whether useful quantities can be computed on 
the lattice to describe the process, which is affected by long-distance contributions that are not 
captured by the Operator Product Expansion (see [jri|] for a more exhaustive discussion). 

From the examples above I conclude that to keep the pace with experiments and help in the 
search of New Physics lattice computations must aim at high precision, typically between a few 
percent and 10% depending on the process and the corresponding non-perturbative parameters 
needed. In order to achieve such an accuracy the computations must start from first principles, 
which implies that the light fermions must be treated as dynamical degrees of freedom, and all 
the systematics associated with renormalization, extrapolation to the continuum limit and chiral 
extrapolation must be kept under control. In the rest of the review I will try to show that each of 
those effects can introduce an uncertainty of 0(5%), and I will do that while presenting a selection 
of recent results for B and D meson decay constants, the Bb(,) parameter and semi-leptonic form 
factors for heavy-light and heavy-heavy transitions. 

In the last part I will describe an approach (not necessarily the only one) in which all these 
systematics can be addressed non-perturbatively. As an application I will present the (quenched) 



computation of the b-quark mass in HQET including 0(l/mb) effects [|T^ 



2. and D^^^ meson decay constants 

In the Standard Model the purely leptonic decays of charged B and D mesons proceed via 
quark annihilation into a W boson. Taking as example the B —>■ TV^ channel, the branching ratio 
can be parameterized as 

^{B- ^T-v,)ocFi\VubW (2.1) 

which turns out to be 0(10^^). The proportionality factor is a function of well-known masses, life- 
times and the Fermi constant. In eq. ( |2.1| ) is the B meson decay constant, which is given by the 
matrix element of the heavy-light axial current between the vacuum and the B-meson state, while 
Vub is the relevant element (actually the smallest and least known) of the CKM matrix. Similarly 

is the non-perturbative matrix element necessary for the SM prediction of the B^ l^^i^ 
branching ratio discussed in the Introduction. 

The B meson decay constant has been computed with three dynamical flavors by the HPQCD 
and the Fermilab, MILC Collaborations ||T3, 14, 15 1. In both cases the rooted staggered quarks con- 



figurations generated by the MILC Collaboration with the AsqTad action have been employed [p^]. 



The heavy b-quark is simulated by using NRQCD in 1 13] and the Fermilab action in [14, 15]. The 



results from |[l3, 14] for = Fe^^/mg^ are shown in figure [I] as a function of the sea quark mass in 
units of the strange quark mass and for the unitary (light sea quark mass equal to the light valence 
quark mass) points only. The curves are the Staggered Chiral Perturbation Theory (S;^PT) JIt] ] 
fits. Although the same formulae have been used and the lattice resolutions are not too coarse, the 
results suggest quite different chiral behaviors (reflected in a 5% difference on the ratio F^JF^), 
probably due to residual cutoff effects. It is interesting to note the consistency of the Fermilab data 
with the curvature predicted from S;^PT, notice however that there the coupling gB*B;r appearing in 
the non-analytic terms has been set to gD*Dn: from the CLEO experiment before performing the fit. 



The final result quoted in 1 13] is Fb = 2 16(9) (19) (4) (6) MeV, where the first error is statis- 



tical (including chiral extrapolations) and the others are estimates of the systematics. The largest 
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Figure 1: SxPT-guided extrapolations for Fq from [O, h4j 



one in particular is due to the matching between the heavy-light current in QCD and in NRQCD. 
This matching involves power divergent mixings between dimension-three and dimension-four op- 
erators in the effective theory and the subtraction has been performed by considering the one-loop 
contribution only. The other systematics included are discretization effects and relativistic correc- 
tions. Most of these cancel in the ratio F^JF^, for which the value 1.20(3)(1) is obtained. 



The Fermilab Collaboration in 014| ] preferred to quote numbers for the ratio only as at that 
time the computation of the relevant renormalization constants was not yet completed. The result 
is FbJFb = 1.27 (2) (6) where the second uncertainty is mainly due to the chiral extrapolation. An 
update including results from two additional lattice resolutions (a = 0.12 and 0.15 fm) and the use 



of the matching renormalization constants computed at one-loop in [ ]18| ] has been presented at this 
conference [|T5|]. The prehminary analysis yields T^b = 191 (5) (8) MeV and FbJFq = 1.30(3) (4), 



both in good agreement with the NQRCD results. 

The ALPHA Collaboration has completed the non-perturbative computation of the renormal- 
ization constant of the static-light axial current with two dynamical flavors in the Schrodinger 



functional (SF) scheme [19]. The main result is the universal (i.e. regularization independent) 
factor <I>(/x)/<I>RGi relating a matrix element ^{jJ.) of the static-light axial current renormalized at 
the scale jj. to its scheme-independent (Renormalization Group Invariant) version. The result is 
shown in figure ^ For /i < 2 GeV perturbation theory fails in reproducing the correct result and 
there would be no way to detect it within perturbation theory only, as the convergence of the series 
appears to be very good in all the range plotted. At the most non-perturbative scale, where large 
volume matrix elements relevant for phenomenology are usually renormalized, the discrepancy 
reaches 5%. The regularization dependent constants needed to match the bare matrix elements 
to the ones renormalized at this scale have also been computed in [^] for different static actions 



(see [ 23 ] for their precise definition) and for the range of bare couplings relevant for simulations in 
large volume using Wilson-Clover fermions. 



4 



SM parameters and heavy quarks on the lattice 



Michele Delia Morte 



*(m)/*rgi 
1.5 



1.4 



1.3 



1.1 



- 


1 ' 




MM| 1 
- 






















SF scheme, 


N,=2 _ 




1 * 






■ 1-loop 7, 2- 


oop /S 




1 


- 2-loop 7, 3- 
1 I I I 


oop jS 

.M.I 



10 100 /. 1000 



Figure 2: Non-perturbatively computed running matrix element of the static-light axial current in the SF 
scheme from []l9[]. The dotted and solid lines are obtained from perturbation theory using 1/2 and 2/3 loops 
expressions for the anomalous dimension of the current [|o| and the j3 -function [pi]]. As an information, the 
A parameter from 1^] is Asf ~ 100 MeV. 



As a first application Fg'^' has been computed on a 24^ x 32 lattice with a ~ 0.08 fm and 
(degenerate) sea quark masses close to the strange quark mass. The result Fg"" = 297(14) MeV 
is rather large compared for example to the quenched value F^^ = 193(6) MeV obtained in p4| ] 
by linearly interpolating in the inverse meson mass between continuum results in the static ap- 
proximation and in the relativistic theory with heavy quarks around the charm. Several effects 
may concur in producing the large Nf = 2 number, for instance cutoff effects, 1/mb corrections 
or sea quark mass effects. While to estimate the latter it is necessary to repeat the computation 
at lighter sea quark masses, for the first two an impression can be gathered by comparing with 
the static result at a similar lattice spacing in the quenched approximation, which turns out to be 
F^\Nf = 0,a ~ 0.08 fm)= 247(5) MeV from @]. This still leaves room for sizeable effects 
of the dynamical fermions, which I will consider again in the following when discussing the Djg) 
meson decay constant. Remaining within the quenched approximation Fb, has also been computed 
including 1/mb corrections explicitly in HQET [p5[]. The final result is nicely consistent with the 
one obtained by the interpolation discussed above, although with larger errors. The computation 
will be described in more detail in the last section. 

Let us now consider the D(s) system. The decay constants Fd and Fd, can be used to extract 
the CKM matrix elements Vcd and Vcs from the CLEO data [^]. The most recent computation 
by the HPQCD Collaboration includes the effects of 2-1-1 dynamical flavors implemented 
in the staggered AsqTad formalism by use of the fourth root of the quark determinant. For the 
valence fermions two different variants of the new Highly Improved Staggered Quark (HISQ) 
action [28] have been used for the light (including strange) and the charm quarks. Some sim- 
ulations parameters are collected in table 1 while results are shown in figure 3, taken from [^]. 
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V 



am,. 



16^ X 48 
20^ X 64 
24^ X 64 
28^ X 96 



0.15 fm 
0.12 fm 
0.12 fm 
0.09 fm 



0.85 
-0.65 
-0.65 
-0.43 



The values of the charm quark mass in units of a are quite 
large in this study and they might cause some concern about 
the size of cutoff effects. In figure 3 however these appear 
to be roughly at the 10% level at the coarsest lattice resolu- 
tion. The concern is then whether all the data are in the scal- 
ing region and a continuum limit extrapolation is justified 
or not. It would therefore be desirable to repeat the com- 
putation at the very fine resolution a = 0.06 fm where the 
MILC Collaboration is indeed producing configurations. 
The final result F^^ = 241(3) MeV, F^JF^) = 1.162(9) is obtained by performing a simultane- 
ous chiral and continuum extrapolation of the data at different quark masses and lattice spacings. 
The overall error includes corrections due to 
the u/d quark mass difference and electro- 



Table 1: Lattice volumes, lattice 
spacings a and values of the charm 

quark mass in units of a from |E7p. 



magnetic effects (see table 2 in [ )27| ] for the de- 
tailed error budget), which make the claimed 
precision clearly impressive. In my opinion 
such a precision calls for a complete clarifi- 
cation of the issues related to the use of the 
"fourth root trick" in dynamical simulations 
of staggered quarks. The discussion on the 
localization, the unitarity and the symmetry 
content of the "rooted" theory [M ^ pll 0] 



> 



is still ongoing and a final conclusion in favor 
or disfavor of it hasn't been reached yet. Also, 
to be able to conclusively judge on the error 
budget it would be useful to have more details 
concerning the Bayesian fits performed, the 
precise functional forms used in the contin- 
uum/chiral extrapolations and also some al- 
gorithmic details. Simulations are indeed de- 
scribed for sea quark masses above one fifth 



0-3 



0.25 



0.25 



0.2 



0.15 




Id 

fa 
U 



0.1 ^ 

0.3 0.4 

Figure 3: Results for the D, Ds (K and n) decay 
constants from for three lattice resolutions 
(see table 1). The chiral fits are performed together 
with those of the corresponding meson masses. 
The continuum limit is given by the dashed lines 
and the final, chirally extrapolated, results are 
represented by the shaded bands. 



of the strange quark mass only [16]. Some of 
these points will probably be clarified in the 
longer publication announced in [^]. 

The European Twisted Mass (ETM) Collaboration has presented at this conference an ap- 
plication in the charm sector of the twisted mass (tmQCD) formalism with two dynamical light 



flavors []33[]. By working at maximal twist the quantities computed are automatically 0(a) im- 



proved [ |34[ ] and no renormalization constants have to be calculated to obtain the decay constants, 
as first pointed out in [|^]. Configurations have been generated for two lattice volumes 24^ x 48 
and 32^ x 64 with lattice spacings a ~ 0.09 and 0.07 fm respectively. The sea quark masses are 
in the range ms/6 and 2m^/3. The decay constants Fd and Fd, have been obtained by interpo- 
lating to the proper value of the meson mass the results produced for heavy quarks around the 
charm. The interpolation in the case of the D meson at the coarser lattice resolution is shown in 
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figure ^ taken from [33]. In tliis case four points liave been fitted with a tfiree-parameters func- 
tional form inspired by HQET. The preliminary results quoted are Fd, = 27 1(6) (4) (5) MeV and 
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Figure 4: Scaling of Fps^/wps as a function of the inverse pseudoscalar meson mass 1/mps. Plot from 



FuJFd = 1.35(4)(1)(7) for a ~ 0.09 fm. The second error comes from the uncertainty on the 
strange quark mass, while the third is due to the uncertainty on the lattice spacing in the case of F^^ 
and to the chiral extrapolation in the case of the ratio FdJFd. The determinations at the finer lattice 
resolution provide consistent results though with larger errors. It is important to assess precisely 
the size of cutoff effects on the result above, as that is obtained interpolating in pseudoscalar meson 
masses which are very close to the cutoff scale (see figure HI). 



Finally, in [36], the QCDSF Collaboration calculated the decay constants of heavy-light pseu- 
doscalar mesons on a very fine lattice (a ~ 0.04 fm, V = 40^ x 80) using non-perturbatively 0(a) 
improved Wilson fermions in the quenched approximation. The result for F^^ is presented in fig- 
ure |5| together with those obtained by the ALPHA collaboration using the same action but in a larger 
range of lattice resolutions [|^, 38 1. The agreement between the results is quite satisfactory and 
suggests the possibility of a joint continuum extrapolation (excluding for example the point at the 
coarsest lattice spacing). The computation of the D meson decay constant requires a chiral extrap- 




Figure 5: The meson decay constant versus a^ with non-perturbatively 0(a) improved Wilson fermions 
in quenched QCD. Results from pq] (star) and [p7[ p^ (squares and circle). 



olation, which in [36] has been performed by linearly extrapolating data corresponding to "pion" 
masses above 500 MeV. An uncertainty associated with this chiral extrapolation is not estimated 
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for the final en^or budget and the values = 220(6) (5) (1 1) MeV (the third error is ascribed to a 
10% ambiguity in the lattice spacing) and fb^/fb = 1-068(18)(20) are eventually obtained. 

The decay constants of B-mesons are also computed in [p6|]. In this case however bare quark 
masses mq with awjq ~ 0.7 need to be considered and the residual, 0{a^), cutoff effects on Fb, are 
estimated by the authors of |36] to be 12%. This sets the limits of the approach. In addition, for 
such masses roundoff effects on the quark propagator at large time separations should be carefully 
checked as well [p^]. 

The different determinations of the D(s) decay constant show statistically significant quench- 
ing effects. For Fo^, which has been computed by most of the collaborations, the results discussed 
are collected in figure |6[ There the errors have been conservatively added linearly. The figure also 
shows the tension, which is emerging with the latest experimental measurement (Fd^ = 275 (10) (5) 
MeV and /Fd = 1 . 24 ( 1 0) (3 ) ) from CLEO-c [ |26| ] . The lattice determinations are indeed system- 
atically below it and in some cases the discrepancy is above two standard deviations. As discussed, 
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Figure 6: Recent determinations of the meson decay constant compared to the experimental result. 



for the system quenching effects appear even larger. This shouldn't be puzzling as for B- 
physics effective theories, rather than relativistic QCD, are simulated. The inclusion of dynamical 
fermions can therefore have different effects in the two cases. 



3. ~^(s) mixing 

The weak interactions induce mixings among flavor eigenstates. At low energies and for B- 
mesons the process is described by the AB = 2 Weak Effective Hamiltonian. In particular the matrix 
elements of four-fermion operators O^^^ (corresponding to the box diagrams) among meson (S^^) 
and S(s)) states need to be computed. The mixing is expressed through the oscillation frequency 

Am(s) 

mBqAmq oc |V,;V,bp(Bq|Ovv+AA|Bq) , (3.1) 

where the proportionality factor is given by the Wilson coefficients (functions of mt/mw and Gp). 
It is customary to introduce the parameter by dividing out the result in the vacuum-saturation 
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approximation 

3 (Bq|OvV+AA|fiq 



Sb„ . (3.2) 



8 F2 ml 

Dq riq 

Oscillations of B mesons are comparatively "slow" and have been observed since UAl, the PDG 
average for Am is 0.507(5) ps^^ On the contrary — mixing is very fast and A/Ms has been 
measured only recently by CDF [pi]], with the result A/Mj = 17.77(10) (7) ps^^ Notice that the 
accuracy of both measurements is at the percent level, which will be very difficult to match from 
the theoretical side. However, by combining these experimental determinations with the lattice 
computations of the Bb(,, parameters the Standard Model values for Vtd and Vts (or ratios thereof) 
could be extracted. 

This year three Collaborations have reported results on the B-parameters with three dynamical 
flavors. The HPQCD Collaboration in [42] has computed B^^ and also the matrix elements for 
AFg on the MILC staggered AsqTad configurations at a ~ 0.12 fm, in a volume 20^ x 64 and for 
sea quark masses equal to one half and one quarter of the strange quark mass. The b-quark is 
treated using NRQCD. The results show very little dependence on the light quark masses within 
the errors and the final estimate isFB^^Bg^^ = 281(21) MeV anAB^^m^) =0.76(11) using two- 
loop formulae for the conversion to the MS scheme. In the computation the operators in QCD are 
related to their NRQCD counterparts including 0(1 /mb) corrections, which bring in operators of 
dimension seven. These operators require a power divergent subtraction, which in 1 42 ] is performed 
at the one-loop level. This means that the subtracted operator is still power divergent. With the one- 
loop value for the coefficient the subtraction itself is about 10% of the final result on Bb., and it gives 
the largest contribution to the systematical error (see table 2 in [p2|]). It is clear that the situation 
becomes worse as finer lattice resolutions are considered, as the subtraction grows linearly with 
I /a. A computation of the subtraction coefficient to higher orders in perturbation theory could at 
least help in reducing the systematic uncertainty associated to the matching. However, as pointed 
out in |p3l], part of this systematic cancels in the ratio t = which can be used to extract W\ 

from . This quantity is now being computed by the HPQCD Collaboration which has presented 
a study using several time sources with smearing to reduce the statistical and fitting errors [43]. 

The Fermilab-MILC Collaborations reported about the work in progress on the computation 
of the ratio E, [ p4| ] employing the Fermilab formalism for heavy quarks and again the MILC con- 
figurations generated at a ~ 0.12 fm. Matching and renormalization (also including 0(1 /mb)) 
are implemented in one-loop perturbation theory. The preliminary results are shown in figure ^ 
(statistical errors only). The light sea quaik mass dependence seems rather small compared to the 
statistical error, whereas the dependence on the light valence quark mass is noticeable within statis- 
tics. To finalize the results the S;^PT formulae for the relevant hadronic matrix elements are being 
determined in order to be able to simultaneously fit the results for different quark masses and lattice 
spacings. Indeed the Collaborations plan to repeat the computation on a finer and a coarser lattice. 

The RBC and UKQCD Collaborations have implemented HQET at the leading order (static 
approximation) combined with light domain wall fermions for a computation of the mixing pa- 
rameters with 2-1-1 dynamical flavors [^]. The lattice used has a linear extent L ~ 2 fm with 
<3 ~ 0.12 fm and Lg = 16, which for the residual mass from the five-dimensional Ward identity 
gives amres = 0.00308(4). Three values of the light sea quark mass have been considered, such that 
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Figure 7: as a function of the valence d quark mass and for different values of the sea quark mass m\. 
Figure from [0- 



the lowest pion mass reached is 400 MeV, while for the static quark the APE and HYP2 [ ]23| ] dis- 
cretizations have been used. The prehminary results F^^^ = 220(32) MeV, F^f/F^''^ = 1.10(;"), 
(nib) = 0.79(4) and Bg^'(mb) = 0.74(10) (in the MS scheme) have been obtained by using 
one-loop mean-field improved estimates of the matching and renormalization factors [ |46| ] and by 
linearly extrapolating the data to the physical point. Large differences between the APE and the 

7 Stat 



HYP2 results have been observed for example for the quantity y ^b^' where the discrepancy 
between the central values is 30%. Notice however that even if a chirally invariant light action is 
used, non-perturbative effects in the renormalization constant of the static-light axial current can 
be large (see figure ^ and in addition static light correlations functions are not automatically 0(a) 
improved, therefore large 0(a) contributions may still affect the results. 

The non-perturbative renormalization programme for the parity-odd static-light four-fermion 
operators in the SF scheme has been completed by the ALPHA Collaboration for the quenched 
case and for two dynamical flavors. In all effective theories the operator Oy™^^ is expanded as 

0^™aaK) =CL(At,mb)Ofv+AA(Ai)+Cs(At,'«b)Oi+Pp(At) + 0(l/mb) , (3.3) 

in other words, already at leading order, and in the continuum, the mixing between the two renor- 
malized operators O'^^+AAil^) ^^'^ ^ss+pp(m) ^'^^ ^'^ considered. On top of that the bare lattice 
operators may mix with operators of the same dimension under renormalization. In particular if 
chiral symmetry is broken by the lattice regularization (like with Wilson fermions) the bare oper- 
ators O^^^^, O'^Y-w^ ^ss+pp ^"^^ ^ss-pp among themselves. In the static approximation, 
it has been shown in [ |47| ] by using symmetry arguments that all the chirality breaking mixings can 
be ruled out if one works with Wilson-tmQCD at maximal twist. ' The renormalization constants 
needed, in a mass independent scheme, can then be obtained by renormalizing the parity-odd op- 



erators C?vA+AV ^sp+sp "^he standard Wilson case where indeed OyA+AV ^'^'^ ^sp+sp ^'^ 
not mix with operators of different chirality. 



' The transformation .'^^■^12 introduced in is not completely well defined. The conclusion is anyway unaffected 
as the absence of mixings in a mass independent scheme can be proven by using the transformations £^'j^i2 and Ex;, only. 
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The non-perturbative universal factor relating a matrix element renormalized at the scale /i 
in the SF scheme to the RGI one is shown in figure || for the operators Oy^'^^y and OyA+AV + 
4C?|pYsp> which renormalize multiplicatively [p9|]. The figure refers to the computation in the 
quenched theory. Perturbation theory seems to work for /X > 1 GeV for these quantities. In par- 



SF scheme, N,= ; 

1-loop 7, 2-loop /! - 

8-loop 7, 3-loop ^ - 
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Figure 8: Non-perturbative running matrix elements for the operators Oy^'^AV (1^^) Oy^^j^p^ ^^^tp+sv 
(right) in the SF scheme (cfr. figure 2). Figure from [p^. 

ticular for c'2 the series might seem badly convergent from the difference between the 1/2 and the 
2/3-loop results, but quite surprisingly the non-perturbative value eventually agrees with the 2/3- 
loop one on all the energy range plotted. Similar findings apply to the Nf = 2 theory as well [pO|]. 
In that case however the final errors on the running of the matrix elements are much larger (up to 
5%), which limits somehow the eventual precision one can reach on the weak matrix elements. An 
improvement might be obtained by repeating the calculation at finer lattice spacings. 




4. Form factors for heavy-light and heavy-heavy semi-leptonic decays 

Semi-leptonic decays of B mesons are still the most precise channel for measuring e.g. |Vub|- 
On the theoretical side they are described in a well-understood way (compared to hadronic decays) 
and experimentally they are easier to study than the less abundant purely leptonic decays. Taking 
as prototype the B nlv transition, the differential decay rate in the SM reads (ignoring the lepton 
mass) 



where q is the lepton pair momentum. The form factor f+{q^) can be extracted from the matrix 
element of the vector current 

{7l{p^)\V^\B{p^))=Mq^){p„+p^+qA^2r+Mq^)q^A,„2 , (4.2) 

with A^2 = (m| — m\)/q^. For vector to pseudoscalar transitions the decay rate is parameterized 
by four form factors which can be obtained from matrix elements of the axial and the vector cur- 
rent 
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The differential decay rate in eq. ( |4.1[ ) grows with the pion momentum and therefore experi- 
mental measurements are more precise for large values of pj^. On the lattice, on the other hand, only 
the low (or large q^) region is safe from large cutoff effects. Incidentally that is also the region 
where HQET is applicable. Notice however that the sensitivity of the matrix element in eq. (12) to 
f+{q^) vanishes for q q^^ax = ("Jb — niji,0) as the kinematical factor in front of f+{q^) vanishes 
in that limit. The form factors are therefore directly computed on the lattice only for some large 
value of q^ and then parameterized over the whole q^ region using functional forms, which include 
kinematical constraints, HQET scaling and dispersion relations as originally proposed in [52]. 

The most recent lattice determination of /+ and /o is due to the HPQCD Collaboration [^]. 
The three-point correlation functions needed to extract the form factors have been computed on 
the same set of //f = 2 + 1 configurations used for measuring Bq^ plus additional sets at lighter sea 
quark masses, down to m\/m^ = 0.125. The b-quark has been simulated in the NRQCD formalism 
with one-loop matching of the currents to 0(l/mb), i.e. including one-loop subtracted dimension- 
four operators. The subtraction in this case contributes about 5% of the final result on the matrix 
element. Four lattice momenta have been used for the pion, p^ = {(0,0,0), (0,0,1), (0,1,1) 
and (1,1,1)} X For each light quark mass the data are interpolated to fixed common values 
of Eji and then extrapolated to the physical point using S;^PT and continuum ;^PT to assess the 
uncertainties in the extrapolation. In the chiral fits the coupling gB*B;r entering the chiral logs 
is left free to vary in order to use the functional form suggested by S^PT also for Eji > Iniji. 
The results are plotted in figure ^ together with the curve obtained from the 4-parameter Ball- 
Zwicky fit [54]. The errors shown are statistical and chiral extrapolation errors only. As expected 



2.5 



• yq-) HPQCD 
■ f^(q") HPQCD 



2 4 6 8 10 12 14 16 18 20 22 24 26 
q' in GeV' 



Figure 9: Form factors (q^) and fo{q^) from [p3[]. Errors are combined statistical and chiral extrapolation 
errors. The curve is the Ball-Zwicky parameterization fit. Courtesy of Junko Shigemitsu. 

from the discussion at the beginning of this section they grow for large q^ and statistic is being 
accumulated to reduce them. The total error in the final budget is 14%, mainly due to statistic, 
chiral extrapolation and matching of the currents. The parameterization of f+{q^) is used to obtain 

/ 2 T^'^'i = 2-0V(41)(39) ps-i , (4.3) 
IVubl ^16GeV2 dq^ 
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which combined with the experimental result from the HFAG [||] for the integrated decay rate in the 
equation above gives |Vub| = 3.55(25)(50) x 10^^. The tension with the inclusive determination 
4.49(33) X 10"^, which in the SM is dis-favored by the global Unitarity Triangle fits Q 



ub 



and poses problems also for Minimal Flavor Violating extensions of the SM []56|], is still there. 

An alternative and complementary approach, which can provide precise results for large q^, 
consists in using Heavy Flavor ;tPT. At leading order in chiral perturbation theory and through the 



order 1 //Mb in the heavy quark expansion [ p7| ] 

2Fn 



*B7r 



1 



(4.4) 



where v is the velocity of the heavy meson (notice [F^*] = 2 in the formula). The method re- 
quires a computation of the coupling gB'B;r> which in the static approximation is obtained from 
the matrix element of the light-light axial current between a B and a B* state at zero momentum 
(and is called g). A very precise determination of g in the quenched approximation has appeared 
this year in [^], while a preliminary A'f = 2 result has been presented at this conference [59]. In 
both cases the HYPl static action p3|] has been used and the required two- and three-point corre- 



lation functions have been evaluated adopting the all-to-all techniques introduced in [60]. In the 
quenched approximation 100 eigenvectors have been computed for the low -mode part of the corre- 
lators whereas in the dynamical case 200 eigenvectors were needed, the number of configurations 
used, on the other hand, was only 32 and 100, respectively. The lattice spacings in both cases were 
quite coarse, 0.1 fm in the quenched computation and 0.2 fm for A'f = 2. The results are collected 



in figure 10 (from [59]) as a function of the pseudoscalar meson mass. In the Nf = case the final 



0.55 

) 

0.5 - 
0.45 ■ 
0.4 



4 Our result (n,-=2) 

O Becirevic et al. (n=2) 

□ Negishi et al. (n^=0) 

A Abada et al. (n =0) 







0.5 



1.5 



mp/(Gev') 

Figure 10: Compilation of results for g = g from [ ^s] , ^ ^ |62| ]. The two error bars on the diamonds 
correspond to statistical and perturbative errors. The empty squares and triangles are non-perturbatively 
renormalized. Figure from 



value g = 0.517(16) is obtained by extrapolating linearly the data in a^m\ while for the prelimi- 
nary dynamical result g = 0.55(1)(3)(3)(6) different extrapolations have been compared. For the 
latter value the first error is statistical, the second from the chiral extrapolation, the third from the 
renormalization factor (computed at one-loop only) and the fourth is an estimate of discretization 
effects. 

Considering now heavy to heavy transitions, the Rome II group in [ ]63[ ] presented a quenched 
computation of the form factor G{w), where w is the scalar product of the velocities of the meson 
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in the initial and in the final state, for the B — > Dlv decay. The computation makes use of the step 
scaling method developed by the group in order to avoid resorting to effective theories. In a small 
volume of linear size Lq = 0.4 fm and a very fine lattice resolution the form factor G{w,Lq) is 
computed at the physical values of the bottom and charm quark masses. This number is of course 
affected by very large finite size effects, which are removed by multiplying (twice) by the step 
scaling function a(L,s,w,mh) defined (for a heavy, would be bottom, quark of mass m^) as 



G{sL,w,mh) 
G{L,w,mh) 



s>l. 



(4.5) 



For L> Lq the step scaling function can not be computed directly at the physical value of the b- 
quark mass, and the key idea of the approach is exactly that it is enough to compute it for ~ ^f^b 
(which typically means that in the last step ~ nic) and then extrapolate it in 1/mh to nib. The 
extrapolation is expected to be smooth as finite size effects (which is what the step scaling function 
describes) shouldn't depend strongly on the heavy-mass scale. This is found to be true in all 
applications of the method (see for recent ones where the step scaling functions have been 
computed also in HQET to turn the extrapolations into interpolations). For the case at hand the 
form factor is finally obtained in a (1.2 fm)^ volume as 



G{w) = a(2Lo, l.5,w,m\,)a{LQ,2,w,rrn,)G(w,LQ 



(4.6) 



Each factor is computed in the continuum limit (although extrapolating in from two lattice 
resolutions only for the a's) and the product is then linearly extrapolated in the light quark mass 
from masses above ms/4. Different values of w have been considered by adopting flavor twisted 



boundary conditions. The result is shown in figure |1 1 




Figure 11: Comparison of |V'cb|G(w) with experimental data [ |65| , p6| . The form factor has been com- 
puted in [M] and the figure has been obtained by extracting |Vcb| at w = 1.2, which yields |Vcb| = 



3.84(9)theo(42)exp X 10-^ Figure from 

The ETM Collaboration also computed the form factors for heavy pseudoscalar to pseu- 



doscalar transitions [67] in the same Nf = 2 setup used for the computation of and Fd, [33] 
(i.e. with heavy quarks around the charm). The preliminary result at a ~ 0.09 fm is consistent with 
the one in []63|] within the still rather large statistical errors. 
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We close this section with the computation of the form factor for B D*lv at zero recoil 
from the Fermilab and MILC Collaborations [^]. This channel is less helicity suppressed than the 
B — > Dlv one and it is therefore preferred for the extraction of |Vcb|. The computation simplifies 
in the zero recoil kinematics as in this limit only one, usually called hp^^ , of the four form factors 
contributes. It is obtained from the matrix element of the heavy-heavy axial current between B and 
D* states. As proposed in [68] the form factor can actually be computed directly from the double 
ratio _ _ 

(D1.y,.|5)(5|j,r,r.y) 3, 

' {D*\cYAc\D*){B\bYAb\B) ' " 

Three similar double ratios had been introduced in [ ]69| ] in order to compute hp^^ (1) to 0(1 /m^) in 
the heavy quark expansion. The expression in eq. (l^/7) gives the correct answer to all orders and 
preserves the feature that most of the lattice current renormalizations cancel in the ratio. Notice 
however that contrary to the double ratios in [|^], has a non-trivial value different from one 
already for mb = rric and therefore the uncertainty on it doesn't strictly scale as — 1- 

In [ ]68| ] the method has been applied on the Nj = 2+1 MILC rooted staggered configurations 
together with the Fermilab formalism for heavy quarks. The results are collected in figure 12. 
The physical, continuum value obtained by using S;^PT formulae for the chiral/continuum ex- 
trapolations is /zai(1) = 0.924(12) (19) where the second error is the sum in quadrature of all 
the systematic ones. By combining it with the experimental measurement (see [Q]), the estimate 
|Vcb| = 3.87(9)theo(V)exp X 10"^ is obtained. 



medium coaise (0.15 fm) 
coarse (0.12 fm) 
fine (0.09 fm) 
extrapolated value 




0.2 

m ' (Gev') 



Figure 12: Form factor /iaj (1) from [ |68| ] as a function of the lightest pseudoscalar meson mass for unitary 
points. The band is the continuum extrapolated curve and the dashed line on the physical point is the total 



error after the inclusion in quadrature of the systematic ones. Figure taken from |68 1. 



5. b-quark mass and B meson decay constant in HQET at 0(l/mb) 

HQET on the lattice was introduced in [^, ^ twenty years ago. It offers a theoretically very 
sound approach to non-perturbative B -physics as it provides the correct asymptotic description of 
QCD correlation functions in the limit m\, —>■ oo. Subleading effects are described by higher dimen- 
sional operators whose coupling constants are formally 0(1 /mb) to the appropriate power. The 
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theory can be treated in a completely non-perturbative way including renormalization and match- 



ing, in principle to an arbitrary order in 1/mb, as it was shown in [J72|]. This implies the existence 
of the continuum limit at any fixed order in the expansion. However precise computations have 
been hampered for a long time by the poor signal to noise ratio in heavy-light correlation functions 
at large time separations, which affects the Eichten-Hill action. The signal can be exponentially 
improved by considering minimal modifications of the action, where the link in the time covariant 



derivative is replaced by a smeared link [73, 23]. The inclusion of dynamical quarks is straightfor- 



ward and the approach can be used together with other methods, as for example the one proposed 



by the Rome II group (see [64] for such recent applications). 



To fix the notation we write the HQET action at 0(l/mb) as 



5hqet = a^Y.i^^'y^^ + ^^)^^ + «spinWi(-aB)vAh + o^kmV^h (-D^) Wi} (5.1) 

X 

with y/h satisfying /+V^h = and P+ = ii^. The parameters Wspin and (0\dn are formally 0(1 /mt). 
For the computation of the b-quark mass the task is to fix 5ya, G>kin and Wspin non-perturbatively by 
performing a matching to QCD. Actually, by considering spin averaged quantities we can imme- 
diately get rid of the contributions proportional to Wspin- I will give here a short overview of the 
computation and present the final results, precise definitions can be found in the corresponding pub- 



lications [ ]12| , [25| ]. Let us start by remarking that in order not to spoil the asymptotic convergence 
of the series the matching must be done non-perturbatively (at least for the leading, static piece) 
as soon as the 1/mb corrections are included. Following [[74|], one can imagine having computed a 
matching coefficient Cmatch for the static theory at order / — 1 in perturbation theory. The truncation 
error ACmatch is 

AC,,atch - [g\rn^)]' ~ I ^, , , 1 » ^ as mb ^ oo , (5.2) 

|^2Zjoln(mb/AQCD) J "Jb 

where is a renormalized coupling at the scale and Bq is the first coefficient of the j8 function. 
In other words the perturbative error due to the matching coefficient of the static term is much larger 
than the power corrections in the large nib limit. In our framework matching and renormalization 
are performed simultaneously and non-perturbatively. 



As the action in eq. ( |5.1[ ) would produce a non-renormalizable theory, we treat the 1/mb cor- 
rections to the static, renormalizable theory as space-time insertions in correlations functions. For 
correlation functions of some multilocal fields O this means 

{O) = (O)stat + 0\ina^Y^{00]dn{x))A3.t + ft'spin«'^ Xl^^'^spin(x))stat , (5.3) 

X X 

where (O)stat denotes the expectation value in the static approximation and Okin(-x) and 0^^m{x) are 
given by \ifh{x)o'R\lfh{x) and i//h(;t:)D^v/h(.x), respectively. We work with Schrodinger functional 
boundary conditions, i.e. we consider QCD with Dirichlet boundary conditions in time and periodic 



boundary conditions in space (up to a phase d for the fermions). For the computation in [12] we 
remain in the quenched approximation. In a small volume of extent Li ~ 0.4 fm, one can afford 
lattice spacings a sufficiently smaller than 1/mb, in such a way that the b-quark propagates correctly 



up to discretization errors of 0{a^). QCD observables defined in this volume are described in HQET 
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2 2 

up to effects of O^^^^^ and O^^^^ . The size L\ is chosen in order to have the two effects 

of the same size. We consider two quantities, <I>^^^(L,mh) defined exploiting the sensitivity of 
SF-correlation functions to the angle B and <I>2™(L,mh), which is given by LT\ where Fi is a 
finite volume effective energy. When expanded in HQET^, <I>f^^'^(L) is given by o\ia times a 
quantity defined in the effective theory (which we call while <I>2*^^^(L) is a function 

of ODkin and mtare = ^m + ^t involving two other HQET quantities (Ff and F^'") . Obviously, 
by equating <I>P*"°(Li,mh) to <I>[^*^^^(Li) one can determine the bare parameters mbare and (O^^ 
as functions of Wh at the lattice spacings used for the volume L\. To eventually compute wib we 
need the phenomenological, large volume, input of the spin-averaged vector-pseudoscalar B-meson 
mass, m^. Here we introduce the step scaling functions Cy{L) to evolve the <I>i's to larger volumes 
and write 

<Q^^(2Li,mh) =£aij(Li)<I>Q™(Li,mh) + 5i2an,(Li) . (5.4) 
j 

Notice that <I>[^*^^^(2Li ,mh) constructed in this way still has a dependence on mh, which is inherited 
from the matching to QCD in L\. The step scaling functions on the other hand are defined in 
HQET and have a continuum limit there. After two evolution steps volumes of extent roughly 
1.5 fm are reached and the bare parameters mbaie and (O^^ can be computed again as functions of 
mh for the corresponding lattice spacings. They are expressed in terms of step scaling functions, 
<I>9*^°(Li,/Mh) and quantities computed in HQET (the large volume version of 0'), F^'^' and 

F^"). At this point the b-quark mass can finally be determined by solving for /Mb the equation 



where = HmL^ooFf andS"^" = -(B|a^i:zOkin(0,z)|S)stat with = 1. In eq. (|3]) 1 have 
emphasized the dependence of the bare parameters mbare and (O^n on mh. However when those are 
re-expressed in terms of <I>P™(Li), ajj, 0',L2), Ff'"(L2) and Fj'"(L2) the equation involves 

only quantities which have a continuum limit either in QCD or HQET. This in particular implies 
that in the procedure all power divergences have been non-perturbatively subtracted. 

To see how the different pieces combine together in the final result it is instructive to consider 
more explicitly the relatively simple case of the computation in the static approximation. In this 
situation only the parameter mbare needs to be determined. In the small volume the matching 
condition reads 

Fi(Li,mh) = Ff (Li) +mbare , (5.6) 

and its large volume version is 

+ '^bare, (5.7) 

to this order we could have just as well used or mB« in the previous equation. If we were able 
to simulate the small and the large volumes at the same lattice spacings we could insert mbare from 



eq. ( p.q ) into eq. ( |5.7| ) and obtain the master equation 

^av ^ (^stat _ Yf\L,)) + Fi (Li ,mh) , (5.8) 



^We set the mass counterterm 8^ in the action to zero here. Its contribution is taken into account in the overall 
energy shift m|,ai-e between the effective theory and QCD. 
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whose solution is the b-quark mass in the static limit. To circumvent the problem we bridge the 
gap in volume in two steps by inserting a step scaling function am (Z. i ) = 2L i (F j''" ( 2L i ) — F j''" (L i ) ) 
into the master equation (i.e. we add and subtract Fj''"(2Li)), which then becomes 

lUm^ - 2Li [E''^' - Ff] - am(Li ) = 2LiFi (Li , m^) . (5.9) 

Now the quantities in the square brackets can be computed at the same values of a and their dif- 
ference has a well-defined continuum limit in HQET because in the combination we are non- 
perturbatively removing all the divergences, particularly the linear one. As announced any refer- 
ence to bare parameters has disappeared in the final equation. 



The graphical solution of eq. ( |5.9[ ) is shown in figure [13|. On the horizontal axis we plot 
: LiMh where Mh is the heavy quark mass in the RGI scheme. The result is M^^^ = 6.806(79) 



□ L,m;''"-L,(E"'-r;"'(Lj))-0_^ 
O 24>,(Lj,M) 



10 



11 



12 



13 



14 



Figure 13: Graphical solution of eq. (5.9) in the quenched approximation. In the caption <I>2(Li,M) 
Liri(Li,Mh). Data from [O]. 



GeV using ro = 0.5 fm to set the scale [75, 76]. 

The inclusion of the subleading 1/mb effects is more involved and I report here the final re- 
sults summarized in table 2 . The different numbers correspond to various matching conditions, 
identified by the choice of the angle(s) 6 and by the strategy adopted, "main strategy" for the first 
line and "alternative strategy" for the second to fourth line. The details can be found in [|^], what 
should be emphasized here is that while there are some differences among the static results de- 
pending on the matching condition chosen, those are completely gone once the 1 / /Mb terms are 
included, signalling practically negligible higher order corrections. The value eventually quoted 
in is mb(mb) = 4.347(48) MeV in the MS scheme. 

As a further application the decay constant of the meson has been computed in quenched 
QCD including 0(l/mb) in HQET [25]. Four quantities are needed for the matching, which again 
has been performed in several different ways. The preliminary results in table 3 show the same 
pattern discussed for the b-quark mass. Notice however that the difference at 0(l/mb) is more 
significant than the errors suggest as most of the uncertainties from the large volume part of the 
computation cancel in the difference. This indeed yields for instance 

(00 = 0, 01 = 1 , 02 = 0) - Fb^^^^^ (00 = 1 , 01 = 0, 02 = 0.5) = 4 ± 2 MeV . (5. 10) 
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00 


roM^"^ 










01 =0 
02 = 1/2 


01 = 1/2 

02 = 1 


01 = 1 

02=0 





17.25(20) 


17.12(22) 


17.12(22) 


17.12(22) 




1/2 
1 


17.05(25) 
17.01(22) 
16.78(28) 


17.25(28) 
17.23(28) 
17.17(32) 


17.23(27) 
17.21(27) 
17.14(30) 


17.24(27) 
17.22(28) 
17.15(30) 


Table 2: Results for the RGI mass Mb from [|l2|]. 


00 


Fl^^^ [MeV] 




+ F^|^[MeV] 






01 =0 
02 = 0.5 


01 =0.5 
02 = 1 


01 = 1 

02 =0 




0.5 
1 


224 ±3 
220 ±3 
209 ±3 


185 ±21 
185±21 
184±21 


186±22 
187 ±22 
185±21 


189 ±22 
189 ±22 
188 ±22 



Table 3: Results for Fb, from [El. 



Finally, the results are in good agreement with the determinations in [|24|,|64|], which also go beyond 
the static approximation. 

6. Conclusions 

A big effort has been devoted in the last years to removing the quenched approximation from 
lattice computations. This is absolutely necessary in order to provide precise theoretical estimates 
to test the Standard Model or to look for signals of New Physics. Several lessons have been learnt 
from these works. Quenching effects have been proven to be large and chiral extrapolations to be 



more delicate than in the Nf = approximation, as partly expected [^, |81[ ]. 

However B-flavor physics is going to become high-precision physics and other systematics 
may significantly affect the results. I have shown that the uncertainties associated to renormal- 
ization, matching, chiral and continuum extrapolations can easily reach the 5 to 10 percent level. 
When choosing an approach for performing a first-principle computation the possibiUty to keep 
these systematics under control should be included among the requirements. 

I have described how these problems can be solved non-perturbatively in Heavy Quark Effec- 
tive Theory on the lattice. The computations I discussed in this framework on the other hand have 
been performed in the quenched approximation only and therefore the results in principle are not 
immediately applicable to phenomenology. The extension to dynamical light fermions is ongoing 
and first steps have been reported at this conference [^]. The method can be used for several 
quantities, the b-quark mass and the B-meson decay constant discussed here but also S — S mixing 
parameters and form factors for semi-leptonic decays. 
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More generally, the non-perturbative matching procedure between HQET and QCD in small 



volume can be adopted also for other effective theories, as it has been done in [|78|, |79|] for a version 
of the Fermilab action. 
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